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Quenched QCD with domain wall fermions 
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We report on simulations of quenched QCD using domain wall fermions, where we focus on basic questions 
about the formalism and its ability to produce expected low energy hadronic physics for light quarks. The work 
reported here is on quenched 8 3 x 32 lattices at fi — 5.7 and 5.85, using values for the length of the fifth dimension 
between 10 and 48. We report results for parameter choices which lead to the desired number of flavors, a study 
of undamped modes in the extra dimension and hadron masses. 



1. INTRODUCTION 

Domain wall fermions [0 were originally pro- 
posed as a technique for putting chiral fermions 
on the lattice. For fermions coupled to a vec- 
tor gauge theory (like QCD), this approach has 
given a way to treat massless (or at least con- 
trollably small mass) fermions, while preserv- 
ing the full continuum symmetry group. Subse- 
quent improvements to the original formulation 
JU made simulations more practical and studies 
of the Schwinger model || , including dynamical 
fermion effects, and quenched QCD g have been 
very encouraging. Here we report a reasonably 
systematic study of low energy, non-anomalous 
QCD physics using domain wall fermions. For a 
review and further references, see j^] . 

2. BASIC QUESTIONS 

For quenched domain wall simulations, there 
are three input parameters for the fermions: to/, 
the explicit four-dimensional bare quark mass; 
L s , the extent of the lattice in the fifth dimension 
and mo, the five-dimensional bare quark mass. 
Except for different symbols for the three param- 
eters listed above, our conventions follow g. 
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A first basic question involves choosing param- 
eters so that the desired number of flavors of light 
fermions is being studied. For free fermions, one 
easily sees that for mo < there are no light sur- 
face states at the ends of the fifth dimension, for 
< to < 2 surface states for a single light quark 
flavor appear, for 2 < too < 4 surface states for 
four light flavors appear and for 4 < too < 6 
there are surface states for six light quark flavors. 
There is a symmetry under too — > 10 — too, which 
remains when coupling to gauge fields is intro- 
duced, although then these values of too shift. 

Figure |l| shows one way to probe the number 
of light fermions. The graph is for quenched sim- 
ulations at (3 = 5.85 on an 8 3 x 32 lattice. For 
a given value of mo, we have measured the chi- 
ral condensate for m/ = 0.02, 0.04, 0.06, 0.08 and 
0.10, and calculated (ifjip(m,f — > 0)) (simple lin- 
ear fits work well). This extrapolation was done 
for three different values of L s , to make sure we 
were close to the large L s limit. The figure shows 
this extrapolated chiral condensate as a function 
of too- 

For mo < 0.9 the extrapolated chiral conden- 
sate is zero, indicating no spontaneous chiral sym- 
metry breaking, due to the absence of surface 
states. For larger values of Too the chiral con- 
densate is non-zero and increasing, until leveling 
off for 1.65 < mo < 2.15. This gives a region 
where there appears to be a single surface state. 
(For free fermions, the normalization of the fields 
depends on mo; this region of insensitivity to Too 
does not appear to have been anticipated.) For 



2 








007 







006 







005 


A 





004 








V 





003 







002 







001 







000 




0.40 1.15 1.90 2.65 3.40 4.15 

m„ 



Figure 1. The four-dimensional chiral condensate 
extrapolated to zero quark mass as a function of 
the five-dimensional domain wall mass. 



larger too the chiral condensate increases as the 
change from a one flavor to a four flavor theory 
occurs. Above about 3.4, a region of relative in- 
sensitivity to too occurs, where there are four fla- 
vors of light fcrmions. The chiral condensate in 
this region is numerically very close to four times 
its value for 1.65 < too < 2.15. 

A second basic question is whether the mixing 
between the chiral modes at s — and s = L s — 1 
can be controlled by choosing L s large. It is 
known that a gauge field with a zero eigenvalue 
for the four-dimensional Wilson Dirac operator 
(with mass too), has no damping for modes prop- 
agating in the fifth dimension. We wanted to ex- 
plicitly measure the effect of this undamped mode 
on non-anomalous low energy QCD physics. To 
this end, we have studied its effect on the chiral 
condensate (for related work see |^|). 

The four-dimensional chiral condensate is de- 
fined through the inverse propagator for domain 
wall fermions, where s = 0(L S — 1) is a source 
for right (left) handed quarks and s = L s — 1(0) 
is a sink for left (right) handed quarks. We have 
studied a generalized condensate, where the dis- 
tance from source to sink is variable. When this 
distance is L s — 1, we have the conventional con- 
densate. 



For a few thermalized configurations, we have 
measured the values of too where there should be 
undamped modes in the fifth dimension. Figure 
|^ shows a measurement of the generalized chi- 
ral condensate as a function of the separation be- 
tween the source and sink in the fifth dimension 
for too very close to the value where a zero of the 
four-dimensional Wilson operator occurs. 
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Figure 2. A generalized version of the four- 
dimensional chiral condensate, where the sepa- 
ration in the fifth dimension between the source 
and sink is s. s = L„ — 1 is the usual condensate. 



One can clearly see the presence of the expected 
translationally invariant mode in the fifth dimen- 
sion. However, this mode contributes about 1% 
to the value of the chiral condensate, for this one 
lattice where its effect should be largest. For an 
ensemble of lattices, at any given too only some 
will have largely undamped modes. (This number 
could be a set of measure zero in the full average, 
although this needs to be studied as a function 
of lattice volume.) Even if undamped modes do 
survive the ensemble average, our example shows 
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that they can make a very small contribution to 
an observable. 

3. HADRON MASSES 

To test the formulation in a more practical way, 
we have measured hadron masses for quenched 
domain wall fermions on 8 3 x 32 lattices with 
L s = 10, 16, 24 and 48 for ft = 5.7 and to = 1.65. 
(Since the conference, we have also simulated at 
L s = 32.) We have chosen m/ in the range 0.02 
to 0.22. For to,/ = 0.02, we have not gotten 
good plateaus with 70-90 configurations and con- 
sequently only report results using to/ > 0.02. 

The rho mass shows no dependence on L s be- 
tween 10 and 48, while the nucleon mass de- 
creases somewhat for larger L s . For L s — 10, 
we find TOp = 0.787(12) + 2.28(5)to/ and m N = 
1.19(3) + 3.71(15)to/. For L s = 48, we find 
TO P = 0.792(28)4-2. 15(12)m/ and m N = 1.12(4)+ 
3.79(19)m/. These are correlated fits to five or 
more valence quark masses and all have \ 2 1 'dof < 
1 . Figure || shows the nucleon and rho masses for 
L s = 48. The to/ = value for for TOAr/m p is 1.42 
(L s — 48), compared to the staggered fermion 
values of w 1.50 0. 
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Figure 3. Rho (□) and nucleon (O) masses versus 
to / for 80 quenched 8 3 x 32 configurations at ft = 
5.7. The fifth dimension has L s = 48. 



We have also measured m T (m/) for these lat- 
tices and find that it is not a simple linear func- 
tion of to/ for the full range of quark masses 
we used. Fits including an m? term, or linear 
fits to the few smallest values for to/, give val- 
ues for to^(to/ — > 0) that agree within errors. 
We then find that mj(0) fits well to the form 
Aexp(-aL s ) + B, with B 0.05. 

We have not yet established an explanation 
for this non-zero intercept for to^.. Our studies 
of the condensate make residual mixing in the 
fifth dimension seem unlikely. Possible causes 
include the finite volume effects which give the 
pion a non-zero intercept in quenched staggered 
spectroscopy, a problem with quenching that is 
brought out by the zero modes of domain wall 
fermions or the observed failure of our integrated 
7T-7T correlator to equal (t/^Vto/ for small values 

of TO/. 

4. CONCLUSIONS 

Our studies of domain wall fermions for 
quenched QCD have shown the technique capable 
of reproducing expected physics. The non-zero 
value for to^(0), even for large L s , deserves more 
attention. Since other formulations have shown 
similar effects for to^(0), we must check again at 
weaker coupling and larger volumes before con- 
cluding that this is a new phenomena revealed by 
domain wall fermions in the quenched approxi- 
mation. 

The authors would like to thank Robert Ed- 
wards, Tim Klassen, Rajamani Narayanan and 
Tom Blum for useful discussions. 

REFERENCES 

1. D. Kaplan, Phys. Lett. B288 (1992) 342. 

2. V. Furman and Y. Shamir, Nucl. Phys. B439 
(1995) 54. 

3. P. Vranas, Phys. Rev. D57 (1998) 1415. 

4. T. Blum and A. Soni, Phys. Rev. D56 (1997) 
174; Phys. Rev. Lett. 79 (1997) 3595. 

5. T. Blum, these pro ceedings. 

6. R. Edwards, et. al. |hep-lat/9802016 . 

7. C. Bernard, et. al. Nucl. Phys. B60A (Proc. 
Suppl.) (1998) 3. 



